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Abstract 

A kinetic theory of relativistic gases in a two-dimensional space is developed in order to obtain the 
equilibrium distribution function and the expressions for the fields of energy per particle, pressure, 
entropy per particle and heat capacities in equilibrium. Furthermore, by using the method of 
Chapman and Enskog for a kinetic model of the Boltzmann equation the non-equilibrium energy- 
momentum tensor and the entropy production rate are determined for a universe described by 
a two-dimensional Robertson- Walker metric. The solutions of the gravitational field equations 
that consider the non-equilibrium energy-momentum tensor - associated with the coefficient of 
bulk viscosity - show that opposed to the four-dimensional case, the cosmic scale factor attains a 
maximum value at a finite time decreasing to a "big crunch" and that there exists a solution of 
the gravitational field equations corresponding to a " false vacuum" . The evolution of the fields of 
pressure, energy density and entropy production rate with the time is also discussed. 
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PACS: 51.10.+y; 98.80.-k; 47.75.+f 



I. INTRODUCTION 

The combination of general relativity with the kinetic theory of gases is remarkably 
useful to construct cosmological models 0. In these formulations the cosmic sources of 
gravitational interactions are represented by the energy-momentum tensor of a fluid; in 
addition we have the hypothesis of homogeneity and isotropy in the form of the well-known 
Robert son- Walker metric ||. Although these theories have explained several important 
features of our universe fundamental questions still remain to be answered ||. 

Models in lower dimensions offer interesting results that, if properly analyzed, can be 
used to gain insight in the realistic formulations. Two-dimensional (2-D) gravity models 
have been under intensive investigation during the last two decades ||, ^ |6|, [7], || [3J . The old 
problem of quantum gravity, black holes physics and string dynamics were tested in these 
theories. In particular Teitelboim Q and Jackiw || proposed a consistent model in two 
dimensions analogous of general relativity. As imediate results, among others 0, this 
model offer a consistent Newtonian limit, gravitational collapse solutions that are basically 
a 2-D Schwarzschild analogue and cosmological models based in a 2-D Robertson- Walker 
metric. 

For cosmological applications, a refinement in the construction of these models can be 
obtained by considering a non-equilibrium scenario, including a bulk viscosity term in the 
energy- momentum tensor (for a review on viscous cosmology up to 1990 one is referred to 
Gron |H|). In the four-dimensional case the inclusion of this term to analyze the evolution 
of the cosmic scale factor with the time was done by Murphy [T1J who has found a solution 



that corresponds only to an expansion. Other models were based on the coupling of the 
Einstein field equations with the balance equations of extended thermodynamics ]T2| (also 
known as causal or second-order thermodynamic theory) and among others we cite the works 
of Belinskii et al JIBf , Zimdahl [T3[] and Di Prisco et al JH| . 



In this work we develop a kinetic theory of relativistic gases in a two-dimensional space. 
The balance laws for the particle flow, energy-momentum tensor and entropy flow are ob- 
tained from the Boltzmann equation. We find also the equilibrium distribution function and 
the expressions for the fields of energy per particle, pressure, entropy per particle, enthalpy 
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per particle and heat capacities in equilibrium in a two-dimensional space. Moreover, by 
using the method of Chapman and Enskog for the kinetic model of the Boltzmann equation 
proposed by Anderson and Witting |[6| we calculate the bulk viscosity and the entropy pro- 
duction rate. We apply the ideas of Murphy [II] to the 2-D gravitational field equations and 
we show that opposed to the four- dimensional case the cosmic scale factor attains a maxi- 
mum value at a finite time decreasing to a "big crunch" and that there exists a solution of 
the gravitational field equations corresponding to a "false vacuum". The difference between 
the solutions in the four- and two-dimensional cases is due to the fact that the relationship 
between the metric tensor and the sources in the 2-D case is modified because the Einstein 
field equations give no dynamics for the 2-D case. 

The manuscript is structured as follows. In Section II we introduce the two-dimensional 
Robert son- Walker metric. The kinetic theory of relativistic gases in 2-D is developed in 
Section III. In Section IV we introduce the gravitational equations of motion in the 2-D case 
and in Section V we search for the solutions of the gravitational field equations. Finally, in 
Section VI we discuss the solutions that came out from the gravitational field equations. 



II. ROBERTSON- WALKER METRIC 



One fundamental feature of 2-D cosmological models is that they show considerably less 
mathematical complexity and at the same time they preserve the physical principles that 
are used to construct their four- dimensional counterparts. One impressive result was that 
in 2-D models the quantization of the gravitational field is consistent || |J, opening the 
possibility of quantum cosmological models for the very early universe. These results in the 
2-D theories are of relevance to include new ingredients in the "realistic" versions |8|. [J . 

As it is well known, the so-called cosmological principle is based on the assumption that 
the universe is spatially homogeneous and isotropic. The metric that describes such kind of 
universe, known as the Robertson- Walker metric, has the following form in a two-dimensional 
Riemannian space characterized by the metric tensor g^ v with signature (H — ) || |9| 

1 — er z 

In the above equation - the so-called cosmic scale factor - is an unknown function of 
the time and has dimension of length, while r is a dimensionless quantity. If we introduce a 
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new variable x = arcsin(y / £r)/ v /£ the equation (TI]) reduces to 

ds 2 = {cdtf - K(tf(dxf. (2) 

The components and the determinant g = — det^g^)) of the metric tensor g^ for the 
Robert son- Walker metric ([!]) with respect to the coordinates (x^) = (ct, x) are 

9oo = 9 00 = 1, 9n = -k 2 = -rr, 9 = (3) 

9 11 

The corresponding non-zero Christoffel symbols read 

r?i = A«, ri = ^. (4) 

where the dot denotes the derivative with respect to the coordinate x° = ct. Once the 
Christoffel symbols are known the non- vanishing components of the Ricci tensor R^ u = RJ lTU 
and the curvature scalar R = g^R^ can be calculated and it follows 

Roo = — , -Rn = —kk, R = 2—. (5) 

K K 
III. KINETIC THEORY OF RELATIVISTIC GASES IN 2-D 
A. Boltzmann equation and fields in equilibrium 

We consider a relativistic ideal gas with particles of rest mass m which is described in 
the phase space by the one-particle distribution function The momentum = 

{p°,p) has a constraint of constant length g^p^p 11 = m 2 c 2 so that f(x^,p^) = f(x,p,t). 

The evolution equation of the one-particle distribution function in the phase space is 
described by the Boltzmann equation, which in the presence of a gravitational field is given 

by mur 



udf_ _ ! df_ _ U^p„ , _ (0) , , 

P dx» P dp C 2 T U 1 h K ' 

where r* is the afline connection. In the right-hand side of the above equation we have 
replaced the collision term of the Boltzmann equation by the model equation proposed by 



Anderson and Witting ||16| , which refers to a relativistic gas in the Landau and Lifshitz 
description [0. Further r is a characteristic time of order of the time between collisions, 
/(°) is the equilibrium distribution function and is the two- velocity such that U^U^ = c 2 . 
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The first moments of the one-particle distribution function are: the particle flow N 11 and 
the energy-momentum tensor T^ v which are defined by 

N" = c[ p"/ V9-, T^ = cf pV7 >/§-. (7) 
J Po J p 

Furthermore, the entropy flow is defined by 

= -he J p»f In fjg^-, (8) 

where k is the Boltzmann constant. 

The balance equations for the particle flow N^, energy-momentum tensor T^ u and entropy 
flux S^ 1 can be obtained from the Boltzmann equation @ and read 

AT^ = , r% = o, S% = <?>0. (9) 

Above q denotes the entropy production rate defined through 

""'^(t -iK JL^s*. (10 ) 



cr J ™ J V/ <0) / / (0) Po 



In the Landau and Lifshitz description the particle flow and the energy-momentum 
tensor are decomposed according to 

= nU* + J», (11) 

= p { ^ u) - (P + w) A"" + ^U»U U , (12) 

where n is the particle number density, the non-equilibrium part of the particle flow, 
pin") jj^Q pressure deviator, i.e., the traceless part of the pressure tensor, P the hydrostatic 
pressure, w the dynamic pressure, i.e., the non-equilibrium part of the trace of the energy- 
momentum tensor, e the internal energy per particle and A Miy the projector defined by 

A"" = g" u - \u»U u , such that A^U U = 0, A"" A" = A" 47 , A^A"" = 1. (13) 
c 2 

The the non-equilibrium part of the particle flow and the pressure deviator p^") are 
perpendicular to £/ M . 

One can also decompose the entropy flow as 

S» = nsU^ + ^ (14) 

where s denotes the entropy per particle and M the entropy flux which is perpendicular to 
If. 
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The equilibrium distribution function in a two-dimensional space, which is the so-called 
Maxwell- Jiittner distribution function, can be written as 

,2 



/ 



(0) 



n 



where 



c 



mc 



(15) 



2mcK 1 (() ' ' * /,V 

where T denotes the absolute temperature. The parameter ( represents the ratio between 
the rest energy of a particle and the thermal energy of the gas. Moreover, K n (() denotes 
the modified Bessel function of second kind (see for example [^UJ ) 



K n (C) 



c 



r(|) 



T(n+i) Ji 



(16) 



Once the equilibrium distribution function is known, one can obtain the following expres- 
sions for the fields of energy per particle, pressure and entropy per particle in equilibrium: 



mc 



^ 2(C) 1 



P = nkT, 



s E = kl--\n 



n 



2mcK 1 (C) 



(17) 



The thermodynamic quantities enthalpy per particle h and heat capacities per particle 
at constant volume c v and at constant pressure c p follows from its definitions by using (|17|), 
yielding 



h = e + 



n 



'dh\ 



mc 



k 



K 2 (Q 



C 2 + 3C 



' de " 



k 



C 2 + 3C 



^2(C) 



' my 

MO. 



c v -\- k. 



' my 



(18) 



(19) 



The above thermodynamic fields in the non-relativistic limiting case where ( ^> 1, i.e., 
for low temperatures, read 



2 kT 
e ~ mc H — — 



1 + ic + - 



h ~ mc + 



3kT 



1 + ic + - 



SE 



T 



— h In V2nmk + In h- + 

2 n C 



+ 2C 4C 2 + 



(20) 



(21) 



In the ultra-relativistic limiting case where ( < 1, i.e., for high temperatures and/or for 
very small rest mass, we have that the leading term of each thermodynamic field is given by 

2kT~ 



P 

e ~kT = — 
n 



Se ~ k 



1 + ln- 



nc 



h « 2kT, 



(22) 
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B. Dynamic Pressure in a Homogeneous and Isotropic Universe 



We shall determine in this section the dynamic pressure and the entropy production rate 
in a spatially homogeneous and isotropic universe. In the four-dimensional case these topics 
were discussed by Weinberg within the framework of a phenomenological theory and by 
Bernstein jlj within the framework of a kinetic theory of gases. Without loss of generality, 
we shall use here the Anderson and Witting model of the Boltzmann equation (ffl) in order 
to simplify the calculations. 

We begin by neglecting the space gradients - since we are dealing with a spatially homo- 
geneous and isotropic universe - and by considering a comoving frame where (U^) = (c, 0) 
in a two-dimensional space. Hence it follows that the Boltzmann equation @ reduces to 

l4- 2 -^ = - 1 (/-/ (0) )- (23) 

OX u K Op CT 

We use now the Chapman and Enskog method (see for example |17|) and search for a 



solution of the Boltzmann equation ( p3|) of the form 

/ = /<°>(l + 0), (24) 

where is the Maxwell- Jiittner distribution function fllpp and f^°'(f) is the deviation from 

equilibrium of the one-particle distribution function. We insert the above representation 

into the Boltzmann equation fl23|) and by taking into account only the derivatives of the 

Maxwell- Jiittner distribution function we get that the deviation <fi is given by 

fdlnfM k d]nf<®\ 
(f) = —CT — — — 2— p- 



dx° k dp 



-CT 



h / K 2 \T c T c kp 2 
n + V~ C K' 1 )f + ¥T VQ T + ¥f~Kp~ Q 



(25) 



For the elimination of h and T from fl25D we use the balance equations of the particle flow and 
energy-momentum tensor of a non-viscous and non-heat-conducting relativistic gas whose 
constitutive equations read: 

N» = nU», = (ne + P)^-- Pg^. (26) 

Insertion of the above representations into the balance equations iV M ;At = and T^ v = 
lead to 
1 d 



(^gnUn = 0, (27) 



«v+ 1 9 



^(ne + P)^^ 



dx v yfg dx u 

In a comoving frame ( p7|) becomes 



+ (ne + P)r 



(28) 



n 



K 
K 



(29) 



while the spatial components of ( p8|) are identically satisfied due to the constraint that all 
quantities n, P, e and k are only functions of the time coordinate. The temporal component 
of (ESI) can be written as 



(30) 



T k k 

where the heat capacity per particle at constant volume c v is given by fllS]^. Equations (E9| 
and ( PD| ) are used to eliminate n and T from (E3I), yielding 



cr 



A; c 
H t^Po 



c p 
kTp 



K 
K 



(31) 



K\ J c v c v T 

Once the non-equilibrium distribution function (p_4j) is known it is possible to calculate 
the projection of the energy-momentum tensor in a comoving frame which corresponds to 
the sum of the hydrostatic pressure with the dynamic pressure, i.e., 



P + zu 



c I p"p u fV9- = cn 3 [ p 2 f 

J p J 



dp 
Po' 



(32) 



We insert fl24j) together with ( pT|) into (|32|) and integrate the resulting equation, yielding 



-7]C- 



where 



K 2 Kii 



1 - 



6 ( . 



(33) 



In the above equation Kix denotes the integral for the modified Bessel functions (see 

p. 483): 

roo g-Ccosht 



Abramowitz and Stegun 



Kin(C) 



-dt. 



(34) 



C J o cosh n t 

Hence we have identified the coefficient of proportionality between w and ck/n as the 
bulk viscosity rj. If we compare (|33|) with the constitutive equation for the dynamic pressure, 
given in terms of the divergence of the two-velocity, i.e., vj = —rjU^.^ we infer that here 
ckjK plays the same role as C/ M !(U . Furthermore due to the fact that the bulk viscosity is a 
positive quantity the dynamic pressure decreases when the universe is expanding [k > 0) 
while it increases when the universe is contracting (k < 0). 
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As in the four-dimensional case the coefficient of bulk viscosity vanishes in the non- 
relativistic and ultra-relativistic limiting cases. This can be seen from Figure 1 where the 



coefficient of bulk viscosity rja/VkTm is plotted versus the parameter ( = mc 2 / (kT). Here 
we have chosen the following expression for the characteristic time |22[ 



1 c 2 c kT 

t = , with v s = J — v —— , (35) 

nav s V c vh 

where a is a differential cross-section and v s the adiabatic sound speed. 

If we consider that the distribution function is given by (^5l), i.e., / = /^(l + 0), we 

can use the approximation ln(l + <fi) « <fi valid for \<p\ <C 1 in order to write the entropy 

production rate fllOD in a comoving frame as 

? = " / f {0) ^V9d P . (36) 



T 

We insert now (|2~3| ) together with ( [31] ) into (|36| ) and get by integrating the resulting equation 

9 , • s 2 

r/c ( k x 



Hence the entropy production rate is connected with the bulk viscosity. Weinberg |2T| has 
derived a similar formula in a four- dimensional space by using a phenomenological theory 
and has also shown that the bulk viscosity alone could not explain the high entropy of the 
present microwave background radiation. For more details one is referred to Weinberg & |21 



IV. GRAVITATIONAL EQUATIONS OF MOTION 

Given the sources of our universe it is possible to obtain the equations of motion for the 
gravitational field in the whole space-time. The starting point is of course, the Einstein 
theory of gravitation. 

As far as we are working in 2-D some fundamental problems appear. In fact, the main 
point is that the Einstein action in 2-D furnishes no dynamics. In other words the usual 
left-hand side of Einstein equations of motion that follows by using the Hamilton variational 
principle 

R^u — -^QpluR = 0, (38) 
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are in fact an identity in 2-D. This is related to the gauge invariances of gravitation in 2- 
D: space-time diffeomorphisms and (local) conformal transformations. With this in mind 
Teitelboim Q and Jackiw |J proposed as 2-D action [^4 



S = J d 2 x^g {N(x) [r(x) + 8ttT^(x)] } 



(39) 



where T^(x) is the trace of the energy-momentum tensor. Above we have not taken into 
account the term that refers to the cosmological constant. Using the variational principle 
for the auxiliary field N(x) the equation of motion that follow is 



R(x) = SttT^x), 

together with the conservation law T^ v = 0. Equation 
with the trace of the energy-momentum tensor. 



(40) 



relates the curvature scalar 



V. SOLUTION OF GRAVITATIONAL FIELD EQUATIONS 

The gravitational field equations are obtained from ( [TOD together with T^ v '. u = and 
by taking into account the constitutive equation for the energy-momentum tensor T^ u = 
eU fJ, U u — (P + w) A^ u where e = ne is the energy density. Hence it follows that 

k = — 47r(e — P — vj s )k ) k(e + P + zu) + ne = 0. (41) 

Since w = —rjk/K, the above system is closed if we can relate the pressure P and the 
coefficient of bulk viscosity r] to the energy density e. Here we follow Murphy Jll| and 
assume a barotropic equation of state for the pressure and a linear relationship between the 
coefficient of bulk viscosity and the energy density, so that 

P = ( 7 - l) e , n = ae. (42) 

Above a is a constant and 7 may range from < 7 < 2. The value 7 = 1 refers to a pure 
matter (dust), 7 = 2 to a pure radiation and 7 < 1 to a "false vacuum". 

If we introduce the Hubble parameter H = k/n the field equations pi| ) can be written 
thanks to (H) as 

H + H 2 = 4vre (7 - 2 -aH), H (7 - aH)e + e = 0. (43) 
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The elimination of the energy density of the two above equations leads to an equation for 
the Hubble parameter that reads 

(H + 2HH) (7 - 2 - aH) + aH(H + if 2 ) + if (7 - 2 - aH) (7 - aH) (if + H 2 ) = 0. (44) 



Let us search for a solution of ( |4"1| ) by considering if = if* =constant. In this case we 
get that (041) reduces to (7 — 2 — ctff*)(7 — aH*)H^ = and apart from the trivial solution 



ff* = we have two other solutions, namely: i) a — 7 /if* and ii) a = (7 — 2) /if*. Only the 
first solution is physically possible since < 7 < 2. 



We proceed to analyze the two energy conditions p3| , namely the weak energy condition 
which dictates that the energy density is a semi-positive quantity, i.e., e > and the strong 
energy condition which imposes that the inequality e + P + w > must hold implying that 
H/H* < 1. 

We are now ready to determine the solutions of the gravitational field equations. For 
that end we write the system of equations (f41~D and (|43| ) as 



(7-l)- = 4vrP( 7 -2-7-), 7- (l --)+£ = 0, (45) 
(7- l)(ij + Jf 2 ) =47rP(7-2-7i7), F 7 (l - H )P + P = 0, (46) 

respectively, where the cosmic scale factor, the Hubble parameter, the pressure and the time 
are now taken as dimensionless quantities with respect to if*. ^From the first system of 



equations (J45|) one may determine the cosmic scale factor K{t) and the pressure P(t) while 
from the second one (|46| ) the Hubble parameter H (t) can be obtained. We have chosen 
two values for 7, one of them 7 = 0.5 correspond to a negative pressure ("false vacuum") 
while the other 7 = 1.5 implies a positive pressure. For the solutions of the two systems of 
equations (f45|) and (|46| ) the following initial conditions were taken into account: 
ffor 7 = 0.5, then if (0) = 0.5, «(0) = 1, P(0) = -1, 

[for 7 = 1.5, then ff(0) = 0.5, k(0) = 1, P(0) = 1. ' ' 

In Figures 2 through 5 it is shown the evolution with respect to the time of the cosmic scale 
factor, the Hubble parameter, the pressure and the energy density which follow from the 



systems of equations (|45"D and (^) by taking into account the initial conditions (flTf ) and the 



barotropic equation of state. 
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Another important quantity which can be plotted versus the time is the entropy produc- 
tion rate. For the cases where 7 > 1 the dimensionless expression for entropy production 
rate ([J7|) can be written as 

< = r^r <«> 

(7 - 1)« 

In order to obtain the above equation we have used the equation of state P = nT and the 
solution of the continuity equation ( |2"D| ) which reads n oc In Figure 6 we show the 

evolution of the entropy production rate with respect to the time. 

VI. DISCUSSION OF THE RESULTS 

We proceed to discuss the results obtained in the last section. First we call attention 
to the fact that the period where this theory can be applied is the one where there exists 
an interaction between radiation and matter. The reason is that in the period of pure 
radiation the dynamic pressure vanishes while in the period of pure matter (dust) there is 
no interaction at all. Furthermore the solutions are not valid for all values of the time, since 
in some period the radiation decouples from matter and we have no more the interaction 
between the radiation and matter which implies a vanishing dynamic pressure. 

In Figure 2 it is plotted the cosmic scale factor as function of the time. We have chosen 
that at the beginning of this period (t = by adjusting clocks) the cosmic scale factor is 
different from zero. We infer from this figure that the cosmic scale factor has a maximum at 
t ~ 0.0243 for 7 = 1.5 and at t ~ 0.0123 for 7 = 0.5. From these points on the cosmic scale 
factor decreases and tends to zero, i.e., it goes to a "big crunch". It is noteworthy to call 



attention to the fact that the corresponding solution in the four-dimensional case [I I] neither 
have this behavior for the cosmic scale factor nor admit a "false vacuum" solution. As was 
previously pointed out the difference between the solutions in the four- and two-dimensional 
cases is due to the fact that the relationship between the metric tensor and the sources in 
the 2-D case is modified because the Einstein field equations give no dynamics for the 2-D 
case. 

The Hubble parameter as a function of the time is shown in Figure 3. For both values 
of 7 the Hubble parameter decrease, attain a zero value for a time which correspond to the 
maximum of the cosmic scale factor and assume negatives values. 
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In Figures 4 and 5 the evolution of the pressure and of the energy density are represented 
as functions of the time, respectively. Both functions decrease with the time and attain their 
minimum values at the times where the cosmic scale factor has its maximum value. Since 
in this theory there exists no mechanism that could increase the pressure and the energy 
density, we infer that the solutions for t^0.0243 when 7 = 1.5 and for t^0.0123 when 7 = 0.5 
are not physically possible. The same conclusion can be drawn out from Figure 6 where the 
evolution of the entropy production rate is plotted as a function of the time. From this figure 
we note that the entropy production rate decreases with the time and attains its minimum 
(q = 0) when the cosmic scale factor reach its maximum value. At this point the entropy 
per particle s attains its maximum value since nDs = q > 0. There exists no mechanism 
in this theory that could increase the entropy density rate from its minimum value with a 
corresponding decrease of the entropy per particle from its maximum value. 
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Figure Captions 

Figure 1: Volume viscosity r/a/^/kTm vs. £. 



Figure 2: Cosmic scale factor vs. time t. 



Figure 3: Hubble parameter H(t) vs. time t. 



Figure 4: Pressure P(t) vs. time t. 



Figure 5: Energy density e(t) vs. time t. 



Figure 6: Entropy production rate q(t) vs. time t. 
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FIG. 3: Hubble parameter H(t) vs. time t. 
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FIG. 5: Energy density e(t) vs. time t. 




FIG. 6: Entropy production rate s(t) vs. time t. 
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